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Abstract. Orthogonal polynomials on the real line always satisfy a three-term re- 
currence relation. The recurrence coefficients determine a tridiagonal semi-infinite 
matrix (Jacobi matrix) which uniquely characterizes the orthogonal polynomials. We 
investigate new orthogonal polynomials by adding to the Jacobi matrix r new rows 
and columns, so that the original Jacobi matrix is shifted downward. The r new 
rows and columns contain 2r new parameters and the newly obtained orthogonal 
polynomials thus correspond to an upward extension of the Jacobi matrix. We give 
an explicit expression of the new orthogonal polynomials in terms of the original 
orthogonal polynomials, their associated polynomials and the 2r new parameters, 
and we give a fourth order differential equation for these new polynomials when the 
original orthogonal polynomials are classical. Furthermore we show how the orthogo- 
nalizing measure for these new orthogonal polynomials can be obtained and work out 
the details for a one-parameter family of Jacobi polynomials for which the associated 
polynomials are again Jacobi polynomials. 



1. Introduction 

The construction of famihes of orthogonal polynomials on the real line from a 
given system of orthogonal polynomials (or from a given orthogonalizing weight fi) 
has been the subject of various investigations [6], [4], [14], [11], [3], [15], [7], [9], 
[16], [25]. Let Pn (n = 0, 1, 2, . . . ) be a sequence of monic orthogonal polynomials 
on the real line, with orthogonality measure /x, then these polynomials satisfy a 
three-term recurrence relation 

(1.1) Pn+l{x) = (X -bn)Pnix) - alPn-lix), n > 0, 

with 6^ G M and > and initial conditions Pq = 1, P-i = 0. The corresponding 
orthonormal polynomials are 

Pn{x) = Pn{x), 

aia2 ...an 
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and they satisfy the three-term recurrence relation 

(1.2) xpn{x) =an+lPn+l{x) + bnPn{x)+anPn-l{x), U > 0. 

Putting the recurrence coefficients in an infinite tridiagonal matrix gives the Jacobi 
matrix 

/bo 1 



J = 



bi 1 



1 



v 



or the symmetric Jacobi matrix 



J.S = 



fbo ai 
ai bi 



a2 



02 

b2 



as 



\ 



Interesting new Jacobi matrices can be obtained by deleting the first r rows and 
columns, and the corresponding orthogonal polynomials are then the associated 
polynomials of order r, denoted by pn\ Instead of deleting rows and columns 
we will add r new rows and columns at the beginning of the Jacobi matrix, thus 
introducing 2r new parameters. We call the corresponding new orthogonal poly- 

( — r) 

nomials anti- associated of order r and denote them by • In Section 4 we will 
give explicit formulas for these anti-associated polynomials in terms of the original 
orthogonal polynomials and the associated polynomials. We also give a fourth or- 
der differential equation in case the original orthogonal polynomials arc classical. 
The construction of the orthogonality measure for the anti-associated polynomi- 
als is done in Section 5. The analysis requires the knowledge of the orthogonality 
measure of the original system and of the associated polynomials, and also the as- 
ymptotic behavior of both P„ and Pn^^ as n — > oo. These things are known for a 
particular one-parameter family of Jacobi polynomials which we will call Grosjean 
polynomials (after Grosjean who gave a number of interesting properties of this 
family in [9]) and which are considered in Section 3. 

The construction of new orthogonal polynomials by changing and shifting re- 
currence coefl&cients has been under investigation by others. If the original system 
of orthogonal polynomials has constant recurrence coefficients, then one essentially 
deals with Chebyshev polynomials of the second kind. Changing a finite number of 
the recurrence coefficients leads to Bernstein-Szego weights [22, §2.6], i.e., weights 
of the form y/l — x'^/p{x), where p is a positive polynomial on the orthogonality 
interval [—1,1]. Such finite perturbations have been considered among others by 
Geronimus [6], Grosjean [7], Sansigre and Valent [21], and Zabelin [25]. The limit- 
ing case, when the number of changed recurrence coefficients tends to infinity (with 
changes becoming smaller) is treated by Geronimo and Case [4], Dombrowski and 
Nevai [3], and in [24]. Instead of starting with constant recurrence coefficients one 
can also start with periodic recurrence coefficients. Finite perturbations of periodic 
recurrence coefficients are considered by Geronimus [6], Grosjean [8] and the limit- 
ing case by Geronimo and Van Assche [5] . Changing a finite number of recurrence 
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Marcellan et al. [11] and the limiting case by Nevai and Van Assche [15]. Associ- 
ated polynomials correspond to deleting rows and columns of the Jacobi matrix, or 
equivalent to a positive shift in the recurrence coefficients. This situation is rather 
well known since it corresponds to numerator polynomials in Pade approximation. 
A general study of such a shift in the recurrence coefficients can be found in Bel- 
mehdi [1] or Van Assche [23], and for Jacobi polynomials we refer to Grosjean [9] 
and Lewanowicz [10]. The case where both a shift in the recurrence coefficients is 
made together with changing a finite number of coefficients is treated by Nevai [14] 
and Peherstorfer [16]. In particular Nevai shows that when the original system is 
orthogonal on an interval A with weight function w, then the orthogonal polyno- 
mials with weight function w{x)/\f+{x)\'^ on A, where f{z) ~ a + S{Bw,z) and 
f+{x) = lime_>o+ /(x + ie), with S{Bw,z) the Stieltjes transform of the function 
Bw and B a polynomial of degree i, have recurrence coefficients which can be ob- 
tained by shifting the original recurrence coefficients and changing a finite number 
of the initial recurrence coefficients. Note that this generalizes the Bernstein-Szego 
polynomials. 



2. Background and notation 



Let us consider the family of monic orthogonal polynomials P^, n = 0,l,2,..., 
where P„ has degree n, defined by the three-term recurrence relation 

(2.1) Pn+i(x) = {x- bn)Pn{x) - a^P^-i(x), n > 0, ^ 0, 

with P_i = and Pq = 1. The sequences 6„ {n > 0) and (n > 1) also 

generate the associated monic polynoimials of order r (r a positive integer) Pn^\ 
n = 0,l,2,...,by the shifted recurrence 



PnU^) = - bn+r)Pi'\x) - al^,p!^\{x), n > 0, 

ir) (r) 

with P_{ — and Pg — 1. The recurrence relation (2.1) can also be written in 
operator form using the non-symmetric Jacobi matrix J given by 

\ 



/ bo 1 

al bi 1 



J 



an 



1 



a: 



r+l C'r+l 



as 

xP = JP, 

where P = (Pq, Pi, P2, •••)*• ^^e same way the recurrence relation defining the 
associated polynomials of order r is represented by 



„Tj(r) _ T(r) fj (r) 
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where P^'^^ = {Plf\ P^'^\ P2^\ ■ ■ - Y and J^^^ is the Jacobi matrix obtained by 

deleting the first r rows and columns of J. 

Observe that J^"^^ — J for every integer r if and only if = & for n > and 
= for n > 0. The corresponding family of orthogonal polynomials are, Tip 

to standardization of the orthogonality interval to [—1,1], the monic Chebyshev 

polynomials of the second kind Un{x) = 2~"C/„(a;), for which b = and = 1/4. 

This fundamental family is a particular case {a = (3 = 1/2) of the Jacobi family 

Pn°''^\x), which are orthogonal polynomials with respect to the weight function 
(1 — x)"(l + x)^ on [—1, 1]. Closely related to the family Un{x) are the Chebyshev 
polynomials of the first kind T„(a:) [a = (3 = —1/2), of the third kind Vn{x) 
{a = 1/2, /3 = -1/2) and of the fourth kind Wn{x) {a = -1/2, /3 = 1/2). These 
monic Chebyshev polynomials satisfy the same recurrence relation 

(2.2) Pn+i{x)=xPn{x)-^Pn-l{x), U > 2, 

but with different initial conditions for Pi and P2- These initial conditions are 
hidden in the Jacobi matrices Jt, Jv, Jw for respectively the polynomials T„, Vn 
and Wni which are explicitly 

where Ju is the Jacobi matrix for the Chebyshev polynomials of the second kind 
Un- These four families can be embedded in the more general two-parameter family 
represented by the Jacobi matrix 

and by definition the corresponding monic polynomials Pn satisfy the recurrence 
relation (2.2) with initial conditions Pi{x) — x — b, P2{x) — xPi{x) — cPo{x). All 
the associated polynomials for r > 1 are the same: 

P^''\x) = Un{x), n>0,r>l. 

3. Grosjean polynomials 

The families {tn-, n > 0} and n > 0} of Chebyshev polynomials have an 
interesting extension given by Grosjean. Let us consider first the monic Jacobi 
family for which the parameters satisfy a + /? = — 1, and define Gq = 1, 

GZ{x) = CnPt'-^-'^Kx), -l<a<0, 

where is a constant making this a monic polynomial, i.e., = 2"^/ {^^~^)- These 
are monic orthogonal polynomials with respect to the weight function 

(3.1) 1 _i<,<i. 

TT \1 + X J 1 + X 

A second family of monic Jacobi polynomials with a + /? = 1 is denoted by 



„o!/'„^ 1 rafcc,!— a) /„\ 



1 ^ „ 
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Withdn = 2''{^''+^) , which are monic orthogonal polynomials with weight func- 
tion 

(3.2) Wgix)^— ^h^— i^ + x), -l<a;<l. 



a{l — a)n \1 + x 

Obviously we have the special cases Gn (x) = tn{x) and gn (x) = Un{x). We 
will refer to the polynomials and g!^ as the Grosjean polynomials of the first and 
second kind respectively. This terminology is justified since Grosjean [9] showed, 
using direct verification, that 

(3.3) (G«)(i) = 



n 1 



and Ronveaux [19] gave the generalization of the differential link t'^{x) = nun-i{x) 
between the Chebyshev polynomials 

(3.4) (G^Yix) = i-ir-\g-^,i-x) = ngltr 

Recall that for a = —1/2 we have Un{—x) = {—l)'^Un{x). Grosjean actually shows 
[9, p. 275] that the only Jacobi polynomials for which the associated polynomials 
are again Jacobi polynomials are the Chebyshev polynomials of the first, second, 
third and fourth kind, and the Grosjean polynomials. Property (3.3) can easily be 
checked from the recurrence coefficients, which are 

(3.5aj ''"-4n2-l' ~ (2n - 1)^ ' 

and 

r, N , -2q! + 1 2 {n + a){n+l-a) „ „ 

(3-5b) ^- = 4(n+ 1)^-1 ' = — — ' 

(see, e.g., Chihara [2, p. 220]) by changing n to n + 1 and a to —a in (3.5a), which 
gives (3.5b). Property (3.4) follows by the differential property {P^'^^)'{x) = 

monic Jacobi polynomials P^'^^ and the symmetry property 

Pn'^\—x) = {—l)'^pil^'°'\x) (see, e.g., Szego [22, p. 63]). It is already interesting to 
note that for the Grosjean polynomials the are rational function of n consisting 
of the ratio of two quadratic polynomials in n (as in the Legendre case), whereas 
in general for the Gegenbauer family one has a ratio of two cubic polynomials in n 
and for the Jacobi polynomials one deals with quartic polynomials in the degree n. 

The pair of Grosjean polynomials also give an answer to the following question. 
Let D denote differentiation D = and let L2 = a{x)D'^ + t{x)D + be the 
second order (hypergeometric or degenerate hypergeometric) differential operator 
for the classical orthogonal polynomials (Jacobi, Laguerre, Hermite, Bessel), and 
{Pn, n > 0} be the corresponding family of orthogonal polynomials, where Pn 
corresponds with the eigenvalue = —n[{n — l)a" + 2t']/2, so that L2{Pn) = 0. 
If L 2 is the formal adjoint of L2, 
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then for which orthogonal family of polynomials {P^, n > 0} does one have 
L2(P*) = 0? Of course, the Legendre polynomials {a' = r) solve this problem 
since then L2 = For Grosjean polynomials the operator L2 is 

LG,a,n = (1 - x^)D^ + (-1 -2a- x)D + n^ L^aAGn) = 0, 

and 

Lg,a,n = (1 - x'')D'' + {l-2a- 3x)D + n{n + 2), Lg,a,n{9n) = 0> 
and we have 

Therefore, if P„ = G^, then P* = g~^i. Prom (3.4) we also have 

D LG,a,n = -f'5,l+a,n-l D- 

The relative position of the zeros x^^^ < X2^n < • • • < x'^^n '^^ compared to 

the zeros xj^n = — cos ^'^^2n^^ -^^ controlled by a classical comparison theorem 
due to Markov [22, Theorem 6.12.2]. The ratio between the Chebyshev weight 

1 1 

wt{x) = - , -l<a;<l, 

TT VI - 

and the weight wq for the Grosjean polynomials given by (3.1) is 

wt\x) I x^ 

— const X ' 



wg{x) \1 — x 

Now (1 + a;)/(l — x) is an increasing function of x on the interval [—1, 1], hence 
for —1/2 < a < the ratio wt/wq is an increasing function on [—1,1], and 
consequently we have the following inequalities 

(2j^l)7r 
2n 

For —1 < a < —1/2 the ratio wt/wq is decreasing, and thus the inequalities for 
the zeros are reversed 

(2j;-^l)7r 
2n 

Similar conclusions can be made for the zeros yj^^ (j ' = 1, 2, . . . , n) of Grosjean poly- 
nomials of the second kind as compared to the zeros — cos (j = 1, 2, . . . , n) 
of the Chebyshev polynomials of the second kind. The ratio of the two weights is 

wu{x) ^ f 1 + x^ 

= const. ' 



x^^ <-cos ^ \^ ' , j = l,2,...,n, -l/2<a<0. 



xl^ > - cos , j = 1, 2, . . . , n, -1< « < -1/2. 



Wg{x) \1 — X 

and thus for 1/2 < a < 2 this ratio is an increasing function so that 

y«„<-cos-^, j = l,2,...,n, l/2<a<2. 
n + 1 

For — l<Q!<l/2 the inequalities are reversed 

yf > - cos j = 1, 2, . . . , n, -1< a < 1/2. 

Finally we note that when we are dealing with Grosjean polynomials of the 
first kind, the product of the weight function of the orthogonal polynomials and 
the weight function of the associated polynomials is constant, as in the case of 
Chebyshev polynomials of the first kind. This is not true when dealing with the 
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4. AnTI- ASSOCIATED ORTHOGONAL POLYNOMIALS 



The situation described in §2 suggests to construct new families of orthogonal 
polynomials, which we will denote by Pn+r , obtained by pushing down a given 
Jacobi matrix and by introducing in the empty upper left corner new coefficients 
6_j (i = r,r — 1, . . . ,1) on the diagonal and new coefficients a?_j ^ {i — r — l,r — 
2, . . . , 0) on the lower subdiagonal. The new Jacobi matrix is then of the form 



' 2 



Ji-r) 



a 



-r+l 



a 



1 

-r+l 
-r+2 



1 

b-r+2 



\ 



an 



\ 



bo 



1 

hi 



J 



We will call the orthogonal polynomials Pn for this Jacobi matrix J^~'^^ anti- 
associated polynomials for the family P^- They contain 2r new parameters and 
satisfy 

[Pn+r i^)]^ — Pn+r (^)- 



For r = 1 we have 



and for r = 2 we have 



h-i 1 
al J 



j(-2) 



a: 



1 
6_i 1 







an 



J 



i — r) 

This new family of anti-associated polynomials PA can easily be represented 
as a combination of the original family P^ and the associated polynomials -P^^i- 
First, denote by Qn the orthogonal polynomials for the finite Jacobi matrix 



^-r-\-l 



1 



V 



a. 



\ 



-1 h-J 



so that they satisfy Qo = l)Q-i=0 and 

Qn+l{x) = (x - h-r+n)Qn{.x) - a^_^j^^Qn-l{x) , 

Then, clearly 



n < r — 1. 



Pt"\x) = Qn{x), 0<n<r. 
For n> r the anti-associated polynomials satisfy the three-term recurrence relation 
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SO that P^~^^J{x) is a solution of the three-term recurrence relation of the original 

family Pn{x). The initial conditions however are Pr~^\x) = Qr{x) and Pj:Zl\x) = 
Qr-i{x), and since every solution of the three-term recurrence relation (2.1) is a 
linear combination of Pn{x) and P^2i{x), we have 

PiZ\^) = ^Pn{x) + BPi'2,{x), n > 0. 

Using the initial conditions for n = and n = 1 gives A = Qr{x) and B = 
—aQQr-i{x), and thus we have 

(4.1) PiZ\^) = Qr{x)Pn{x) - alQr-i{x)P!;^\{x), n > 0. 

From this representation it is easy to construct a differential equation satisfied 
by the family P^^ ^^e original family P„ is itself solution of a linear differen- 
tial equation of second order, for instance when are the classical polynomials 
(Jacobi, Laguerre, Hermite, Bessel), then they are a solution of the hypergeometric 
differential equation 

(4.2) L^y = a{x)y" + T{x)y' +(^{n- l)a" - r'n^ 2/ = 0. 

The techniques used in [17], [18], [19], [20] and the fact that 

(4.3) L*PW^ = {a" - 2r')P; 

easily give a fourth order differential equation satisfied by the anti-associated poly- 
nomials P^ZJ {x). The general technique is the following. Let B{x) = —aQQr-i{x) 
and put J{x) = B{x)P^2^{x). Then we can transform the equation (4.3) to 



L*,^-{a"-2r')P;„ 
and introducing the differential operator 
i?2 = (tB^D^ + [{2a' - t)B^ - 2aBB']D 

+ 2a{B'f - aBB" - BB'{2a' - r) - [^(n^ - n - 2) + t'(1 + n)]B^ 

this is equivalent to 

(4.4) i?2 J = {<j" - 2t')B^P'^. 

From (4.1) we then have 

R^P^-;} = R2\QrPn\ - - 2T')alQl_^P^, 
and eliminating the second derivative using (4.2) then leads to 



D T>{~^) T^/T T> I AT r>> 
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where Mq and A^o are polynomials. Taking the derivative in (4.5) and using (4.2) 
to eliminate also gives 

where Mi and A^i are polynomials, and repeating this also gives 



{-r)ili/ 



n+r 



M^Pn + N2P 



This shows that 



(4.6) 



det 



(-r) 



\a[a[R2Pt+M M, nJ 

which is the desired fourth order differential equation. When Pn{x) = G^{x), then 
a" — 2t' = so that (4.4) simplifies and becomes homogeneous. The differential 
equation (4.6) however remains one of the fourth order, except when a = —1/2, 
because then 

Pi-l}{x) = Ao{x)Tr,{x) + Boix^x), 

where Aq and Bq are polynomials. Similar as in the above reasoning we then get a 
second order differential equation 

/ PiTl^ Ao Bo\ 

Ai Si =0, 

A2 B2J 

which is equivalent to the equation given in [21]. 



Mo No\ 

Ml Ni 



= 0, 



det 



p(-r) 
n-\-r 

\cT[cT[pi7Jyy 



5. Construction of the orthogonality measure 

If we use probability measures throughout the analysis, and if we use lower case 
p and q for the orthonormal polynomial, then we have pn = ^nPn, where 

In = {aia2 ■ • •an)~^ 

Similarly 

Qn{x) 



qn{x) = 



n < r. 



and Pnli = In^Pnli where = ai^n- Thus the orthonormal anti-associated 
polynomials are 

pi^ri^) = {a-r+i ■ ■ ■ ao)"Sn^'n+? (a;), 
and using this in (4.1) gives 



(5.1) 



pi+r{x) = qr{x)pn{x) -qr-l{x)p^^li{x) 



ai 



The orthonormal polynomials are useful in obtaining the weight function for the 
anti- associated polynomials. Indeed, we can compute the weight function using 
Christoffel functions 



K{x) 



p]{x) 
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Theorem MNT. Suppose Pn{x) are orthonormal polynomials for a measure fi on 
[—1,1] and let jj. belong to Szego's class, i.e., 



Then 

(5.2) lim nXn{x) = 7Tfi'{x)\^l — 

n — >oo 

holds almost everywhere on [—1,1]. 

The Szego condition can be relaxed and in fact it suffices to assume that is a 
regular measure on [—1, 1] (i.e., supp(//) = [—1,1] and lim^^oo Tn^" = 2) and 



/ 



b 

log iJ,'{x) dx > — oo 



in order that (5.2) holds almost everywhere on [a, 6] C ( — 1, 1) [12, Thm. 8]. We 
will assume these conditions for the measure /i and moreover we allow the addition 
of a finite number of mass points to fx. Then (5.2) will still hold almost everywhere 
on [a, b] . Indeed, if we add a mass point c to fx then by the extremum property 



An(x;//)= min ql{t)dn{t), 

qn{x)=l J 



where the minimum is taken over all polynomials of degree at most n which 
take the value 1 at the point x, we see that for the measure //c = A* + ^^c (i-e., the 
measure fx to which we add a mass point at c with mass e) 



K{x;i^c)^ min [ j ql{t)dii{t) + eql{c) 
> min / ql{t)dn{t) 



= Xn{x;ii), 

so that liminf„^oo 'T'An(a;; /i) > 7r//'(a;) Vl — x"^, almost everywhere on [a, 6]. On 
the other hand, consider the polynomial qn{t) = {t — c)rn-i (t) /{x — c), where rn-i 
is the minimizing polynomial for the measure di'{t) = (t — c)^ dii{t), then for a; 7^ c 

K{x;i^c) < 7 — ^-T2 I cfrl_i{^)dii{t) 
[x - cy J 

= {x-c)-^ j rl_,{t) dvit) 
= (x- c)~'^Xn-i{x;i') 

so that limsup^^;^ nXn{x; jjc) < li^n^oo{x — c)~^nA^(a;; v) = 7r/i'(a;)\/l — x"^ al- 
most everywhere on [a, 6], since u is regular and satisfies the Szego condition on 
[a,b]. Combined with the previous inequality this gives 

lim nX^(x: — ttij,'(x)\/'\ — x"^ 
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almost everywhere on [a,b], independent of the mass point c. This procedure can 
be repeated, so that adding a finite number of mass points does not change the 
behavior in (5.2). 

In the previous discussion we wanted to allow the measure fj. to have a larger 
support than [—1, 1] (by allowing mass points outside [—1, 1]) but in such a way 
that (5.2) still holds. Alternative we can use a weaker result by Nevai [13, Thm. 54, 
p. 104] in which the support of /j, is allowed to be [—1, 1] U E, where E contains at 
most a denumerable number of points which can only accumulate at ±1. 

Theorem N. Suppose n G M(0, 1), i.e., the recurrence coefficients an and bn have 
asymptotic behavior given by 

lim Qn = 1/2, lim b^ = 0. 

Then 

limsupnAn(a:; /i) = 7r/i'(a;)\/l — x'^ 

n— »cx3 

holds for almost every x e supp^n) . 

With this result, if we can compute the limit of nA^(x; ^) for x G (—1,1), then 
this limit is almost everywhere equal to 7r|U'(x) yl— x^, which allows us to compute 
the weight /x' on (— 1, 1) without having to worry about the fact that supp(//) may 
be larger than [—1,1]. The existence of such a limit is however not guaranteed 
by Nevai's theorem, contrary to the theorem of Mate, Nevai and Totik where the 
existence of the limit is in the conclusion of the theorem but where supp(//) = [—1,1] 
is required. 

If we assume that the original family pn belongs to a measure /x in the class 
M(0, 1), then the anti-associated polynomials p^^"^ will also have a measure 
which belongs to the class M(0, 1), and thus we know that has support [—1, 1]U 
E, where E is at most denumerable with the only accumulation points at ±1 [24, 
Thm. 1, p. 437]. In fact, there can be at most 2r mass points, since the associated 
polynomials of order r of the anti-associated polynomials pn are again pn and 
these have all their zeros inside [—1, 1], and the interlacing property of orthogonal 
polynomials and associated orthogonal polynomials shows that adding one row and 
one column in the Jacobi matrix to form the anti-associated polynomials pn can 
at most add a mass point to the left of —1 and to the right of 1. Adding r rows and 
columns thus can add at most r mass points to the left of —1 and r mass points 
to the right of 1. If the original measure belongs to the Szego class (which is a 
subclass of M(0, 1)), then the measure /i^"^) restricted to [—1, 1] also belongs to the 
Szego class. This is so because if we denote by A^fZ^^j the restriction of to the 

interval [—1,1], then lim„^oo 7„(/x*^~^))/7„((Uj_[''-^j) exists and is strictly positive 

[13, Thm. 25 on p. 136] and since 7n(A*^~'^'*) = 7n(A*)(<^-r+i " " '0'q)~^ and n belongs 
to the Szego class, which is equivalent with the statement that lim^^oo exists 
and is strictly positive, it follows that lim^^oo 7n(A*[r[\])/2"' exists and is strictly 

positive, so that lJi)^_-^ ^ belongs to the Szego class. 

In order to use the result in (5.2) we observe that for n > r 

n r— 1 n—r 
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and using (5.1) this gives 

n—r n—r / \ 2 n—r—1 

3=Q k=0 ^ fe=0 

- — gr(a:)gr-i(2^) E^'fe(^)^'fe-i(^)- 

In order to be able to compute these sums, we will need to be able to compute 
Christoffel functions of associated polynomials and sums of mixed form containing 
the product pn{x)p^^}_i{x) . This is in general not so easy. However, when we 
use Grosjean polynomials of the first kind, then the associated polynomials are 
Grosjean polynomials of the second kind, so that we are always dealing with Jacobi 
polynomials, for which we can compute these sums, at least as n becomes large. 
For the anti-associated polynomials corresponding to Grosjean polynomials we thus 
can prove the following result. 

Theorem 1. If the original system of orthogonal polynomials consists of Grosjean 
polynomials of the first kind, then the anti-associated polynomials pn are orthog- 
onal with respect to a measure which is absolutely continuous on [—1, 1] with 
density 



Wr{x) 



sin(— Tra) (1 — x)"" 
TT (l + a;)"+i 



qr{x) - aoe''"''qr-i{x) 



{1-xY 



(l + x) 



a+l 



-2 



In addition, there may be at most 2r mass point outside [—1, 1] (r to the left of —1 
and r to the right of 1 ). These mass points are the roots of the equation 

\x — ll"" 

qr{x) - aoqr-i{x) sign x . i^+i = 0- 

Proof. If the original family consists of Grosjean polynomials, then by (4.3) and 
(3.1) we have for — 1 < x < 1 

nl^ I = VI - x^ sin(-.a) ^\l-lr ' 

and this even holds uniformly on closed subintervals of ( — 1, 1), and similarly by 
using (3.2) (with a replaced by —a) 

A\ 1- 1 '^V^\ (1)/ m2 1 -2a(l + a) (l-x)" 

(5.4) hm - > \pi'(x)r= , 7 ---^ ---r- 

n^oon ^ Vl -a;2 sin(-7ra) (l+x)«+i 

For the remaining sum, which consists of a mixture of the original polynomials and 
the associated polynomials, we can use Darboux's extension of the Laplace-Heine 
formula for Legendre polynomials [22, Thm. 8.21.8 on p. 196] 

-a-l/2 . n\ -/3-1/2 

V^Pt'^^ (cos 9) = (^sin - J (^cos - J 

(5.5) COS f fn + (a + /3 + l)/2l^ - ^^±i/^7r^ + 0(l/n), 



UPWARD EXTENSION OF A JACOBI MATRIX 



13 



which holds uniformly for x = cos 6 on closed intervals of ( — 1,1). For Grosjean 
polynomials Pn{x) = G'^{x) we have 



Pn{x) = [l + 0{l/n) 



2n7r 



sin(— ttq;) 



-mra{l + a ) j-,(-a,i+a) 
sin(— ttq;) 



n-l 



and hence for x = cos 9 with e < 6 < tt — e 

-a-l/2 



Pnix) = 



2 9 

— sm - 

sm(— ttq;) \ 2 



cos 



a+1/2 



C^+1/2 
COS nf) TT 



+ 0(l/n) 



and 



Pn-l(a^) = 



sin(— ttq;) \ 2 



sm 



cos 



COS ( nO + 



a-l/2 



-TT 



+ 0(l/n). 



Using this gives 

r 1 ^ \ (1) ^ ^ V-2a(l + a) 

hm - > pfc a; p^^i a;) = ^ 



1 



sin(— TTCt) sin 6*72 COS 6*72 



X 



1 / 

lim — y COS I 



Using 2 cos a cos h — cos (a + 6) + cos (a — h) gives 

i g,o, (u - ^-^.\ cos - ! 'I\ 



k=l 



^ n—r 

— [cos (2k9 - 7r/2) + cos an] 
2n ^ 



k=l 



and thus uniformly for x on closed subsets of (—1,1) we have 

1 "S? ^ \ (1) ^ \ \/-2a(l + ")cos7ra 
(5-6) hm - ypk{x)pl^^{x) = — . n ■ r \ • 

n^oo n Smfc'sml — TTQ!) 



fe=l 



Alternatively, (5.5) can also be obtained by using a Christoffel-Darboux type for- 
mula [1, corollary 2.12] which for orthonormal polynomials is 
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When the original system consists of Grosjean polynomials of the first kind one 
knows by (3.3) that the associated polynomials are Grosjean polynomials of the 
second kind g~" and by (3.4) we know that the derivative of is ngl^t-i- For the 
orthonormal polynomials this gives 



Pn{x) = [l + 0{l/n)]. 

pi'\x)^[l + 0{l/n) 

p'^{x) = [l + 0{l/n)]. 

and hence using (5.5) gives 

1 " 

lim -y2Pk{x)Pkli{x) 

n—>-oo n ' 



2n7T 



sin(— ttq;) 





i + o) 


sin(— 


7ra) 


/ 27m 









sin(— ttq:) 2 



Pt'-'-'^\x), 
P^-"'i+«)(x), 

(l+a,-a)/ X 
1-1 y-^'Jj 



k=l 



-2a(l + Q;) , 

lim 



sin(— TTCt) sin n^oo 



. c^ + 3/2 \ / -a + 3/2 
cos [n + \)d — — TT 1 cos I nd — - — tt 



a + 3/2 \ / -a + 1/2 
— cos I nv — — TT 1 COS [n + 1)6^ — - — tt 



Simple trigonometry then gives (5.6). Combining the limiting relations (5.3), (5.4), 
and (5.6) gives uniformly for x on closed intervals of (—1, 1) 



1 " 1 



1 



J=0 



Vl - x'^ sin(-7rQ;) 

ao 



ql{x) 



(1 + a;) 



a+l 



(1-x)' 



— 2 — gr(a:^)?r-i(a:^) cos Q;7r-\/— 2q;(1 + a) 



Observe now that a\ = 2q:(1 + a), thus we find 



lim - Yj^^r^\x)f 



n— >oo 77, 



j=0 



1 (1 + X 



sin(-7ra) ^/T^^ (1 - x)^ 



qr{x) - aoe''°"'qr-i{x) 



(l-x) 



(l + x) 



a + l 



and this holds uniformly on closed intervals [a, b] C (—1, 1). From (5.2) we can then 
conclude that the orthogonality measure for the anti-associated polynomials 

has an absolutely continuous part on ( — 1, 1) with weight function 



(5.7) Wr(x) = 



sin(— ttq;) (1 — xY 



Qrix) - aoe"^''qr-i(x)——, — — r 
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which is the weight function of the original system Pnix) (Grosjean polynomi- 
als of the first kind) divided by a positive factor containing the new parameters 
b-r, . . . , b-i and o?_^_|_i, . . . , ai^, Oq. This factor cannot vanish on (—1, 1), because 
then both the real part and the imaginary part of 



need to vanish. The imaginary part can only vanish when qr-i{x) — 0, and assum- 
ing this, the real part can then only vanish when also qr{x) = 0. This is impossible 
since two consecutive orthogonal polynomials cannot have a common zero. 

Since the original orthogonal polynomials are Jacobi polynomials, they will be- 
long to the class M(0, 1). Moreover, the original system satisfies 

oo 

I](|l-4al+i| + 2|6fc|)<oo, 

fe=0 

and hence also the new system of anti-associated polynomials satisfies this trace 
class condition. But then we know [24, Theorem 6] that the orthogonality measure 
IJ,^~'^^ is absolutely continuous on (—1, 1) and wc have obtained the weight function 
in (5.7). The mass points outside (—1, 1) are those points a; G M \ (—1, 1) for which 

J2'k=o\Pk '^^]^ ^ This means that at a mass point x we have Pn{x) — > 0, and 
from (5.1) this implies that 



lim Pn{x) 



Qrix) qr-l{x) 

ai pn{x) 



For X ^ [—1, 1] we know that |Pn(^)| increases exponentially fast, hence at a mass 
point we always have 

Qrix) qr-i{x) lim — — = 0. 

ai n^oo pn{x) 

The limit of the ratio Pn-ii^) / Pn{x) can be found by using Markov's theorem, 
which states that this limit is the Stieltjes transform of the measure /i, or it can be 
obtained from Darboux's generalization of the Laplace-Heine formula for Legendre 
polynomials for x outside [—1, 1] [22, Theorem 8.21.7]. Both methods give 

.(1) 



lim ^^-1^"^^ - 



n^oo aiPn{x) {x + 1)"+^ 



where the right hand side is to be taken positive if a; > 1 and negative if a; < — 1. 
This means that a mass point satisfies 

(x — 1)" 

Qrix) - aoqr-i{x)j^^-^-j^ = 0. 
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Remark: In order to determine the measure //*^~'") for the anti-associated poly- 
nomials corresponding to Grosjean polynomials, one can also use the results from 
Theorem 3.9 in Pcherstorfer [16] using the Cauchy principal value and the Stieltjes 
transform. Indeed, the Grosjean polynomials are one of the seldom cases where 
a nice explicit expression for the Cauchy principal value and the Stieltjes trans- 
form exists and for this reason the results of [16] can be applied without prob- 
lems. Observe that this technique is basically also the one used by Grosjean in [9] . 
Our approach using Christoffel functions has the advantage that it advoids taking 
boundary values of a Stieltjes transform or evaluating a Cauchy principal values 
and uses only information of the orthogonal polynomials on the interval [—1, 1]. In 
particular our method also works when appropriate asymptotic information of the 
orthogonal polynomials and the associated orthogonal polynomials on the real line 
is available. 



6. Examples 

The simplest examples occur when we take a = —1/2, in which case the original 
system consists of Chebyshev polynomials of the first kind. The weight function 
for the anti-associated polynomials then becomes 



,,11 1 1 Vl 

Wr{x) — 



TT y/l-x^ \qr{x) - iaoQr-i/ Vl - x^\'^ TT (1 - x^)q^{x) + alq^_^{x) ' 

Hence this weight function is the weight function of Chebyshev polynomials of the 
second kind, divided by a polynomial which is positive on [—1, 1]. Such orthogonal 
polynomials are known as Bernstein-Szego polynomials [22, §2.6]. 

In case the original system consists of Chebyshev polynomials of the second 
kind, for which all the recurrence coeflBcients are constant, we have pn = Un and 
p^n-i — Un-1- Relation (5.1) then becomes 

pl^+ri^) = qr{x)Un{x) - 2aQqr-l{x)Un-l{x) . 

Using Un{x) = sin(n -|- 1)9/ sinO {x = cos^), one easily shows for —1 < a; < 1 



1 1 

lim -Yuf(x) = 

n^oo ' 2(1 - 



J(l-x2) 



and 



and hence 



1 

n. — >rv-i n ' * 



X 



n^oon^ ' ^ ' 2(1 -a;2)' 



1 " 1 

- ^(a^)]^ = 2(1 -x'^) ^^^^^^ ~ 4.aQxqr{x)qr--i{x) + 4ag?^_i(a;)] 



J(l -a;2) 

and hence the weight function becomes 



2 ^fl 



X 



2 
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i.e., this is again a Bernstein-Szego weight. Observe that it can be written as 

W{X) = ■ 2- 

\qr{x) - 2aoe'-^qr-i{x)\ 
For the mass points we see that they can occur only for x ^ [—1, 1] when 

qr{x) - 2aQqr-i{x) lim ^'^-i(^) = q, 

n—>-oo Un{X) 

and since Un-i{x)/Un{x) — >• l/{x + \J x'^ — 1) for x ^ [—1, 1], where the limit is to 
be taken positive for x > 1 and negative for a; < —1, it follows that x is a mass 
point only when 

(6.1) qrix) - 2aQqr-x{x) ] = 0. 

a; + va; — 1 

If we multiply the left hand side by qr{x) — 2aoqr-i{x) / {x — y/x^ — 1), then this 
implies 

q'^{x) - 4aoxqr{x)qr-i{x) +4aog^_i(a;) = 0. 

So the mass points are zeros of the polynomial in the denominator of the weight 
function, but only those zeros for which (6.1) holds. 
For r = 1 we can consider the following cases: 

(1) ao = 1/^2 and 6_i = 0. In this case the weight function of the anti- 
associated polynomials is w{x) — l/(7rVl — x'^) and thus we have the 
Chebyshev polynomials of the first kind. 

(2) ao = 1/2 and b-i = —1/2. Then the weight function is the one for Jacobi 
polynomials with a = 1/2 and /3 = ~l/2 and thus we have Chebyshev 
polynomials of the third kind Vn{x). 

(3) Finally when ao = 1/2 and 6_i = 1/2 the weight function becomes the 
one for Jacobi polynomials with a = —1/2 and P = 1/2, so that we have 
Chebyshev polynomials of the fourth kind Wn{x). 
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